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Abstract 

We analyze the quantum Bianchi I model in the setting of the nonstandard loop quantum 
cosmology. Elementary observables are used to quantize the volume operator. The spectrum of the 
volume operator is bounded from below and discrete. The discreteness may imply a foamy structure 
of spacetime at semiclassical level. The results are described in terms of a free parameter specifying 
loop geometry to be determined in astro-cosmo observations. An evolution of the quantum model 
is generated by the so-called true Hamiltonian, which enables an introduction of a time parameter 
valued in the set of all real numbers. 
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I. INTRODUCTION 



The great challenge is quantization of the Belinskii-Khalatnikov-Lifshitz (BKL) theory 
[l|-[3j. It is a generic solution of General Relativity (GR) that does not rely on any symmetry 
conditions. It presents an evolution of the universe near the space-like cosmological singu- 
larity (CS) with diverging gravitational and matter fields invariants. The singular solution 
may be applied both to the future singularity (Big Crunch) and past singularity (Big Bang). 
The BKL scheme also appears in the low energy limit of superstring models, where it is 
linked to the hyperbolic Kac- Moody algebras [4|. 

The dynamics of the BKL model, close to the singularity, may be approximated by the 
two key ingredients: (i) Kasner-type evolutions well approximated by the Bianchi I model, 
and (ii) spike-type transitions described by the Bianchi II model. The former case occurs 
when time derivatives in the equations of motion are important, whereas the latter case 
appears when space-like gradients play the crucial role. The general BKL type evolution, 
near the cosmological singularity, consists of a sequence of epochs (i) and (ii), which may 
lead to the oscillatory and chaotic type dynamics. 

Present paper concerns quantization of the Bianchi I model with massless scalar field. It 
is a companion paper to our recent paper \^ , presenting classical dynamics of the Bianchi I 
model in terms of loop geometry. 

Through the paper we apply the reduced phase space, RPS, quantization method devel- 
oped by us recently (6|-10]. It is an alterative method to the Dirac quantization recently 



applied to the quantization of the Bianchi I model |llHl6 



In section II, we recall some elements of the classical formalism |5[ for self-consistency. 
We redefine an evolution parameter and elementary observables, and introduce the true 
Hamiltonian. Section III is devoted to the quantization of the classical model. Compound 
observables are quantized in terms of elementary observables. We present solution to the 
eigenvalue problem for the volume operators. Examination of an evolution of the quantum 
system completes this section. We conclude in the last section. Unitarily non-equivalent 
representations of the volume operators are briefly discussed in the appendix A. We make 
comments on the operators ordering problem in the appendix B. 



II. PREPARATIONS TO QUANTIZATION 

The Bianchi I model with massless scalar field is described by the line element 



ds 2 = —N 2 dt 2 + £ a 2 (t) dx 2 , (1) 



i=i 

where 

k 



-) , dr = Ndt, £k 1 = l = J2 k * + k l> ( 2 ) 



and where k^ describes matter field density (k^ = corresponds to the Kasner model). For 
a clear exposition of the classical singularity aspects of the Bianchi I model, in terms of the 
loop geometry, we recommend 12| . 
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A. Hamiltonian constraint 



The gravitational part of the classical Hamiltonian, for the Bianchi I model with massless 
scalar field, reads pf 



H g = -7" 2 / d 3 x Ne- l e l3k E^E bk F l ab , (3) 
Jv 

where 7 is the Barbero-Immirzi parameter, V C S is the fiducial volume, £ is spacelike 
hypersurface, N denotes the lapse function, e^k is the alternating tensor, Ef is a densitized 
vector field, e := ^\ det E\, and where F* b is the curvature of an SU(2) connection A l a . 

The resolution of the singularity, obtained within LQC, is based on rewriting the curvature 
F k b in terms of the holonomy around a loop by making use of the formula [5[ 

where 

h Uij = h^hfi\h^)-\h^ ) y x (5) 

is the holonomy of the gravitational connection around the square loop D^-. The loop is 
taken over a face of an elementary cell, each of whose sides has length fijLj with respect 
to the flat fiducial metric °q a b '■= $ij °u l a °u)l; the fiducial triad °e\ and cotriad °u k satisfy 



°uj % „ "e - = 8j] Ar Dij denotes the area of the square D^-; LiL 2 L 3 = Vq, where Vq = f v y/°qd 3 x 



a 



is the fiducial volume of V with respect to the fiducial metric. 

The holonomy in the fundamental, j = 1/2, representation of SU(2) reads 

h { t ] = cos(^Ci/2) I + 2 sin(/i iCi /2) n, (6) 

where Tj = —ioijl (a - , are the Pauli spin matrices). The connection A k a and the density 
weighted triad E% are determined by the conjugate variables Ck and pu as follows 

A i a = a L- 1 X, El = Pi Lf 1 L," 1 X, (7) 

where 

c i = ^d i L i) \pi\ = dj a k Lj L k , (8) 
( the dot over a, denotes derivative with respect to the cosmological time), and where 

{Ci,pj} = 8nG~f5 i:j (9) 

Making use of (J3J) and (j4]) leads to H g in the form [|3| 

F„ = lim #W^3) ; ( 10 ) 

where 

(11) 

and where V = aia^asVo is the volume of the elementary cell V with respect to the physical 
metric q ab = a x a 2 a z °q ab . 
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The total Hamiltonian for Bianchi I universe with a massless scalar field, 0, reads 



H = H g + « 0, (12) 

where H g is defined by ({TO]) . The Hamiltonian of the scalar field reads = 
^pl|PiP2P3| _5 /2, where and p^ are the elementary variables satisfying {(f), p^} = 1. The 
relation H ~ defines the constraint on phase space of considered gravitational system. 

Making use of fl6]) we calculate (jTTj) and get the modified total Hamiltonian corre- 
sponding to (1121) in the form 



H w /N = 1 sgn(pip 2 P3) 

87TG7 2 /il/i 2 /i3 



sin(ci/ii) sin(c 2 /i 2 ) /i 3 sgn(p 3 )^/ -|— - =- + cyclic 



2V 

(13) 

where V = \Z\piP~2P3~\- In what follows we assume that 




^k-=\r^K (14) 



where A is a free parameter of our model. 

The choice (fl4"|) for leads, in the Dirac quantization [13|, |l4j, |l6|, to the dependance of 
dynamics on V C E. In what follows, we specialize our considerations to the case V = S = 
T 3 . In such a case the volume observable, V — f v ^d 3 x (where g denotes the determinant 
of the physical metric q a b on £), characterizes the entire space part of the universe. Thus, 
V is chosen unambiguously and V is physical. In the case when one considers the Bianchi I 
model with the M. 3 topology, the volume V C S = M 3 is only an auxiliary tool devoid of any 
physical meaning 1 . 

The present paper is a quantum version of our recent paper [5}, where we consider the 
Bianchi I model with the T 3 topology. The aim of both our papers is presenting a quantum 
Bianchi I model in terms of the nonstandard LQC, which is an alternative to the standard 



LQC results [14j, |l5[ (with V = T 3 and the choice (JHJ). The results obtained within these 
two methods are similar. Detailed comparison is beyond the scope of the present paper, but 
will be presented elsewhere after we complete quantization of the Bianchi II model [18| |. 

Let us analyze the dependance of our results on the choice of coordinates in S, and 
consequently on the choice of the fiducial volume Vq = LiL 2 L 3 . It is clear that is a 
coordinate length, whereas a^L^ is the physical one. The latter is invariant with respect to 
a change of the system of coordinates so we have a^Lk = a' k L' k , while Lk — » L' k . Since, due 
to (jHJ), we have 

Cj = 7— aiLi, \pi\ = ajLj a k L k , (15) 

the canonical variables q and Pi do not depend on the choice of the coordinate system. 
Thus, the destination variables, defined by Eq. (120]) . share this property too. The holonomy 



It is tempting to introduce the notion of a local volume, i.e. a sort of 'density' of the physical volume 
operator V. It appears that the local volume may be defined in the case of any topology of space, however 
it would depend on the choice of coordinates on E. Similar notion has been used by Martin Bojowald 
(see, Appendix in [l3|)> while considering lattice refining in loop quantum cosmology. 
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variable h k ^ depends on /ifcCfc, and with our choice (TH|) leads to 



c k 



A' 



Pk\ 



(16) 



which proves that the holonomy variable does not depend on the choice of coordinates 2 . The 
flux variable is pk so it does not depend on the choice of coordinates either. Since holonomy 
and flux are basic variables of the formalism, our final results do not depend on the choice 
of coordinates (the choice of Vq). This is why the variables (3k and Vk share this property as 
well. Our results do depend on V, but it is correct since V = T 3 is the whole space. 

We wish to emphasize that f fl3|) is not an effective Hamiltonian for quantum dynamics 
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but a classical Hamiltonian modified by approximating the curvature F* b by holonomy 
of connection around a loop with finite length. Our approach is quite different from the 
so-called polymerization method where the replacement c — > sin(c/i) / p, in the Hamiltonian 
is treated as some kind of an effective quantization. Our method has been presented with all 
details and compared with the Dirac quantization method in 0, Q, M, 13] • For an extended 
motivation of our approach we recommend an appendix of jsj. 

In the gauge N = V = A/IP1P2P3I5 the Hamiltonian modified by loop geometry reads 



sgn(pip 2 )|piP2| 3/2 sin 




sin (A 2 ) + cyclic 



+ f- (17) 



Since we consider the relative dynamics [!, 0] our results, in what follows, are gauge inde- 
pendent. 

Equation (jl7jl corresponds to the effective quantization of the standard LQC In the 
nonstandard LQC Eq. (j!7p is treated as the constraint, which is to be imposed into the 
classical dynamics. In the standard LQC one implements this constraint into an operator 
constraint defining quantum dynamics (kernel of this operator is used to find the physical 
Hilbert space). Thus, in the reduced phase space quantization (nonstandard LQC) there is 
no quantum Hamiltonian constraint, contrary to the Dirac quantization, i.e. the standard 
LQC (which is motivated from LQG). In both cases one modifies, to some extent, gravity 
theory: already at the classical level in the reduced phase space quantization, only at the 
quantum level in the Dirac quantization. 

The Poisson bracket is defined to be 



{•,-}:=8vrG 7 



where (ci, c 2 , c 3 ,pi,p 2 ,P3, • 
phase space is defined by 



fc=i 



d- d- 



d- d- 



dc k dp k dp k dc k 



d- d- 



d- d- 



d(p dp^ dp^ dcj) ' 



(18) 



are canonical variables. The dynamics of a function £ on a 



i {^H (x) }, £ G {c 1 ,c 2 ,c 3 ,p 1 ,p 2 ,p 3 ,(j),p 4> }. 



(19) 



The dynamics is defined by the solutions to f fT9|) satisfying the constraint w 0. The 
solutions of ffl9l) ignoring the constraint are nonphysical. 



We have A = A' since A is a physical length. 
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In what follows we shift from (ck,Pk) to another canonical variables (vk,Pk) 

Pk := -^=, v k := |p fc | 3 / 2 , (20) 

VlPkl 

(where k = 1, 2, 3) since they are proper variables to examine the singularity aspects of our 
model In this paper we restrict our considerations to Vk > 0, since we wish to ascribe 
to it a directional volume observable. 



B. Correspondence with FRW observables 

We redefine the original elementary Bianchi observables that has been found in (BJ as 
follows 

^v^m (21) 

and 

1 /|tan(^)|\ 3 sgn^O. + Ofc) 



1 3k V \ / 2^^0102 + 0x03 + 0203' V ' 

where k 2 := 47rG/3. Dropping subscripts leads to the elementary observables of the FRW 
model found in [9|. One may verify that the algebra of redefined observables is 

{O 4 ,O,} = 0, {^,0,} = ^-, {A,,A J } = 0. (23) 

Our main concern is quantization of the volume observable defined as follows 0] 

V = M 2 ^ 3 ) 1/3 , (24) 

where 

, , / 3VttG (Oj + O k ) <f> \ , . 

vt = 3k 7 A \OA cosh v 3 JY - 3kA< . 25 

In an 'isotropic' case (i=j=k) we get the expression for the volume observable of the FRW 
model 



C. Redefinitions of evolution parameter 

Since the observables Oi are constants of motion in and <f> £ R, it is possible to make 
the following redefinition of an evolution parameter 

(p := (26) 

2^0102 + 0x03 + 0203 

so we have 

Vi = 3«7A|Oi| cosh 3k ((Oj + fc ) y? - A*), (27) 
which simplifies further considerations. 
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D. Redefinitions of elementary observables 

One can make the following redefinitions 

A := Ai - (Oj + k ) if. (28) 
Thus, the directional volume (127]) becomes 

Vi := \wi\, Wi = 3firyAOj cosh(3K*4.j). (29) 
The algebra of observables reads 

{0;, 0,1 = 0, {A i ,O j } = 5 ij , {Ai,Aj} = 0, (30) 
where the Poisson bracket is defined to be 

{ '' ' } := 2 \dA~kdO~k ~ dO k dAk)' (31) 
fe=i 



E. Structure of phase space 

All considerations carried out in the previous section have been done under the assump- 
tion that the observables Oi, O2 and O3 have no restrictions. The inspection of ( I2"2"j) . ( 125]) 
and (12 8j) shows that the domain of definition of the elementary observables reads 

D:={(A k ,O k )\A k eR, O 1 O 2 + O 1 O 3 + O 2 O 3 >0}, (32) 

where k = 1,2,3. The restriction Oi0 2 + O1O3 + 2 3 > is a consequence of the 
Hamiltonian constraint (see, [5| for more details). 
In what follows we consider two cases: 

1. Kasner-unlike dynamics: (a) Oj > 0, Oj > 0, O k > 0, which describes all three 
directions expanding (b) Oj < 0, 0-, < 0, 0^ < 0, with all directions shrinking. 

2. Kasner-like dynamics: (a) Oj > 0, Oj > 0, Ofc < 0, which describes two directions 
expanding and one direction shrinking; (b) Oj < 0, Oj < 0, 0^ > 0, with two directions 
shrinking and one expanding. 

This classification presents all possible nontrivial cases. Our terminology fits the one used 
in [l3[ due to the relation Oj = 6nkiK, (0 < K = const), where constants ki are defined 
by®. 



F. True Hamiltonian 

Now, we define a generator of an evolution called a true Hamiltonian H. Making use of 
(128|) . and Oj = const (see we get 

{Ai, H} := ^ = -(Oj + O k ), {O t , H} := ^ = 0. (33) 

dip dip 

The solution to (1331 is easily found to be 

H = 1 2 + Oi0 3 + 2 3 . (34) 
The true Hamiltonian is defined on the reduced phase space which is devoid of constraints. 
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III. QUANTIZATION 



A. Representation of elementary observables 

We use the Schrodinger representation for the algebra (1311 defined as 

~ h d - 

fe -)• O k f k (x k ) := - - — /fc(x fc ), A -)■ A k fk{xk) ■= x k fk(x k ), k = 1,2,3. (35) 

One may verify that 

[64, 6,1 = 0, [A,^] = 0, [A, 6,]=^, (36) 

The representation is defined formally on some dense subspaces of a Hilbert space to be 
specified later. 



B. Kasner-unlike case 

The condition O1O2 + O1O3 + O2O3 > is automatically satisfied in this case, because 
Oi,02 and O3 are of the same sign. To be specific, let us consider (la); the case (lb) can 
be done by analogy. 

Let us quantize the directional volumes by means of Wi defined in (1291) . A standard 
procedure gives 3 



w :- 



(6 cosh (3kA) + cosh (3kA) 6 J = - — (2 cosh(fcr)— + bsmh(bx)) , (37) 
2 V J 2 \ dx / 



where a := K'yXh and b := 3k, and where we have used the representation for the elementary 
observables defined by (I3l)|) . 

In what follows we solve the eigenvalue problem for the operator w and identify its domain 
of self-adjointness. 

Let us consider the invertible mapping L 2 (M, dx) 3 ip — > Uip =: / G L 2 (I, dy) defined by 

,/,n n ItP-V&f^n 

Uj,{z):= ^ \ /2 ^l ,U =:f(y), ^R, yeI:=(0,7r/6). (38) 
sin 1 (by) 



We have 







= / ipip dx 
J —00 

|tg 1/6 (^) |tg^(^) In |tg^(^) |)) 
?(ln|tg 1 / fe (^)|)^(ln|tg^(^)|). 



2 JUTy 10 y 2 J, sm(by) 



l *H**(*)l)«M**(*>l)^ = <& ^ (39) 



/ sin 1 / 2 ^) sin 1 / 2 ^) 
In what follows we drop subscripts of observables to simplify notation 
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Thus, the mapping (138|) is isometric and hence unitary. 

Now, let us see how the operator w transforms under the unitary map (138]) . The trans- 
formation consists of the change of an independent variable 

2 

x i — y y '.= -arctan(e 6:r ), (40) 

which leads to 

ia ( d , . . ,, .\ . d .ab . , , . 

zcoshfox)- — hosmh(ox) h- > — ia- — h i — cot(oy), (41) 

2 V dx J dy 2 

and re-scaling with respect to a dependent variable 

d 



— la- 
dy 



+ i— cot(by) i — y sin 1 ^ 2 {by) ( — ia- — h i— cot(oy) J sin 1 ^ 2 (6y) = —ia—. (42) 
z y ay Z J ay 



In the process of mapping 

w w U' 1 = —ia— =: w, (43) 
dy 

we have used two identities: sin(fry) = 1/ cosh(6a;) and sinh(6x) = — cot(by). 

Since w > (for O > 0), we assume that the spectrum of w consists of positive eigen- 
values. To implement this assumption, we define w := yftiP and consider the eigenvalue 
problem 

~ a2 ^ fu = u2fu ' y E ^ n ^- ( 44 ) 

There are two independent solutions for each value of v 2 (where v £ R), namely: sin(^y) 
and cos(^y). Removing this degeneracy leads to required positive eigenvalues of w. We 
achieve that in a standard way by requiring that the eigenvectors vanish at the boundaries, 
i.e, at y = and y = ir/b. As the result we get the following spectrum 

f u = Nsm(-y), v 2 = (nab) 2 , n = 0,1,2,... (45) 
a 

It should be noted that for n = 0, the eigenvector is a null state and thus the lowest 
eigenvalue is v 2 = (ab) 2 . Next, we define the Hilbert space to be the closure of the span of 
the eigenvectors f T45]) . The operator w 2 = — o 2 ^^is essentially self-adjoint on this span by 
the construction. Due to the spectral theorem [19j we may define an essentially self-adjoint 

operator w = y— & 2 Jpr as follows 

wf u := uf u , v = ab, 2ab, 3ab, . . . (46) 

We have considered the case w > 0. The case w < does not require changing of the 
Hilbert space. The replacement w —w leads to v \-> —v. 

Finally, we find that the inverse mapping from L 2 (I, dy) to L 2 (R, dx) for the eigenvectors 
of w yields 

>v s , s ~ i , x , s sin ( 2 xarctg(e 6x )) , 

sin -y) = /„(y ^ CT 1 /* y := ^(x) = U *\ ^ - (47) 
a cosh y/ (fa) 
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C. Kasner-like case 



In the case (2a), the conditions Oi0 2 + O1O3 + 2 3 > with Oi < 0, 2 > 0, 3 > 
are satisfied in the following domains 4 for 

O x G(-di,0), 2 G(d 2 ,oo), 3 G(d 3 ,oo), (48) 

where d 2 > di, and where d 3 = d\d 2 j[d 2 — d\) so cf 3 > d\. The full phase space sector of the 
Kasner-like evolution is defined as the union 

|J (-di,0) x (d 2 ,oo) x (d 3 ,oo) (49) 

Q<d 1 <d 2 

In the case of 2 and 3 , the restrictions for domains (|48|) translate into the restrictions 
for the corresponding domains for the observables w 2 and w 3 , due to (1291) . and read 

w 2 e(D 2 ,oo), w 3 e(D 3 ,oo), (50) 

where D 2 = n^\d 2 and D$ = KjXd^. Thus, quantization of the w 2 and observables can 
be done by analogy to the Kasner-unlike case. The spectra of the operators w 2 and 11)3 are 
almost the same as the spectrum defined by (|46l) with the only difference that now v > D 2 
and v > D 3 , respectively 5 . 

The case of Wi requires special treatment. Let us redefine the elementary observables 
corresponding to the 1-st direction as follows 

ocosh(Mi) 

The transformation (151 p is canonical, since {Qi,Q 2 } = 1, and invertible. The domains 
transform as follows 

OiG(-di,0), AiER — > fii G (0,di/6) =: (0,Di), fi 2 G R. (52) 
The observable v\ in terms of redefined observables reads 

v x = < ^n 1 (i + s%), ^6(0,00), (53) 

where ab/H = ^nG'jX. To quantize observables O x and Q 2 we use the Schrodinger repre- 
sentation 

^2 -> Ckfix) := -ihd x f(x), fix -> h x f(x) := xf(x), f G L 2 (0, D x ). (54) 

Let us find an explicit form for the operator y (Qi + corresponding to (1531) . Since 

Qi > 0, the following classical equality holds 

flrflj = Q k -Q 2 - Q{- k - m ■ Q 2 ■ fij 1 , (55) 



4 The case (2b) can be done by analogy. 

5 Spectra are insensitive to unitary transformations. 
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where m, k £ R. This may lead to many operator orderings at the quantum level. This issue 
is further discussed in the appendix B and in the conclusion section. 
We propose the following mapping (we set h = 1) 

Q1Q2 -> SH| := ~ ^*fi 2 ^^^ m fi 2 ^ + ^fi 2 ^^"^2 ^ = -xS^ - d x + rnkx- 1 , 

(56) 

which formally ensures the symmetricity of The second equality in (I56I) may be 

verified via direct calculations. 

Now, we define the following unitary transformation W 

L 2 ([0,D l ],dx) 3 f(x) 1 y Wf{x) := e L 2 ([0,2^/D;],dy). (57) 

One may verify that we have 

wd x w^ = -d y - 1, ly^^t = 4%, - 4^ + A ■ (58) 
y v v v y 

Thus, the operator W transforms (156]) into 

- 8 ty + ^( 4mk -\)- ( 59 ) 
The eigenvalue problem for + reads 

-C + ^( 4 ^-i)+x) $ = ^- (6°) 

Now, we can see an advantage of the chosen ordering prescription (|56|) . It enables finding a 
very simple form of the volume operator. Taking k = m = 1/4 turns (160]) into 



2/ 



2 



4 



The problem is mathematically equivalent to the one dimensional harmonic oscillator in a 
'box' with an edge equal to 2\fD\. There are two independent solutions for a given v 



^^-^(-i.+ii^, 

4 4§,£). 



(62) 
(63) 



where i-F\ is a hypergeometric confluent function, $^1 and $^2 are even and odd cylindrical 
functions, respectively. A standard condition for the symmetricity of the operator defining 
the eigenvalue problem ( 16TI) leads to the vanishing of the wave functions at the boundaries 
(as the box defines the entire size of the 1-st direction). The solution (after retrieving of h 
and ab) reads 6 . 

*=^-M-^4lfD- (64) 



We ignore the solution <£>„,i because it cannot vanish at y = 0. 
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The solution ( 164]) vanishes at y = as $ is an odd function. The requirement of vanishing 
at y = 2y/Di leads to the equation 

flu 3 3 2DA . . 



An explicit form of (1651) reads 



E ( " m +!), f^y = °. < 66 > 



(I) V ft , 

n=0 v2/n v ' 

where (a) n = a(a + 1) . . . (a + n — 1). It results from (|6"6"]) that the eigenvalues must satisfy 
the condition: z/ > ab. 



D. Quantum volume operator 

Classically we have 

V = \w lW2 w 3 \ 1/3 . (67) 
One may verify that Poisson commute and commute, so we can take 

V 3 := V1V2V3 = \wiW2Ws\. (68) 
The eigenfunctions of the operator W1W2W3 have the form ir A i. A 2,A3 ._ 
/i Al ( x i)/2 A2 (a ; 2)/3 A3 (2 ; 3), where fi Xi (xi) is an eigenvector of u>j with eigenvalue Aj. 
The closure of the span of _F Al ' Aa ' A 3 is a Hilbert space, in which V 3 is a self-adjoint operator 
(by construction). 

Due to the spectral theorem on self-adjoint operators [l9|, we have 

V = (V 3 ) 1/3 > yp^U^2,\a . = □^ A 1, A 2, A 3 ; (69) 

where 

□ := |AiA 2 A 3 | 1/3 . (70) 



1. Kasner-unlike case 

In the Kasner-unlike case we use the formula (146]) to get 

D = \riin 2 n 3 \ 1/3 ab, n x , n 2 , n 3 G 1, 2, 3, . . . , (71) 

which shows that the spectrum of the volume operator does not have equally distant levels. 
The volume □ equal to zero is not in the spectrum. There exist a quantum of the volume 
which equals A := ab = 12irG , -fXh, and which defines the lowest value in the spectrum. 



2. Kasner-like case 

The spectrum in this case reads 

□ := [J D dl>d2 , D dl:d2 := {\ dl \ d2 \ d3 \d 3 = d 1 d 2 /(d 2 - d^}, (72) 

0<di<d 2 

where A^ is any value subject to the condition ( 166|) . X d2 > D 2 and \ d . A > D 3 are given by 
( 146]) . The volume □ equal to zero is not in the spectrum. 
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E. Evolution 



In this section we ignore the restrictions concerning the domains of Oi, O2 and O3, and 
we assume that the Hilbert space of the system is L 2 (M. 3 ,dxdydz). An inclusion of the 
restrictions would complicate the calculations without bringing any qualitative change into 
the picture of evolution. 

The generator of evolution determined in (1511) may be formally quantized, due to ([Ho*]) , 
as follows 

H^H= -h\dyd z + d x d z + d x d y ). (73) 

Since it is self-adjoint in L 2 (IR 3 , dxdydz), a quantum evolution can be defined by an unitary 
operator 

JJ = e -itir{d y d z +d x d z +d x dy) ^ r G R. (74) 

Let us study an evolution of the expectation value of the directional volume v\ 

(Vltrtit/lV) (75) 
Since V\ does not depend on y and z, we simplify our considerations by taking 

jj^ _ e -ihT(d z +d y )d x _ (JQ^ 

If we are interested in the action of U\ on the functions f(x) G L 2 (IR, dx), then the derivatives 
— i-^y- and — i-^ occurring in U\ commute and, being self-adjoint, lead finally to real numbers. 
Let us call them k y and k z , respectively, and let us introduce the parameter k = k y + k z . 
Hence, U\ further simplifies and reads 

U l =e khTdx . (77) 

The action of U\ on fix) reads 

U 1 f(x) = f(x + khr). (78) 

We recall that under the unitary mapping L 2 (R,dx) ^ L 2 (I,dy), defined by flM}, the 
operator vi becomes —io,-^ on L 2 (I,dy). Now, let us study an action of operator Ui on the 

functions ip(y) G L 2 (I,dy). Straightforward calculation leads to 

i, , , s v(l arctan(e 6x )) o , , , , 

L 2 (I, y) 3 tpiy) ^ V^L) e L 2 (R, x), (79) 

and we have 

arctan(e te )) arctan(e fe3:+fe ^ r )) 

^i— T7^ = wo 1 (8°) 

cosh ' (bx) cosh ' (bx + bkhr) 

The transformation U~ l gives 

<^(|arctan(e ba;+6Mr )) <^(| arctan(e 6fc/lT tan(^))) 

— — 775 ^ — : =: <f T {y), (81) 

cosh (6x + fiAjftr) sin (fy)(tan(^)e 6fe ^ + cot(%)e- bkhT ) 

where ip T= o(y) = <p(y). Now, we observe that the symmetricity condition 

(<p T (y)\vi<p T (y)) = (vivMlvAy)) (82) 
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leads to 

^r(fK(f)-^r(0)Vr(0) = 0. (83) 

We use the result (IHTI) to calculate the limits 

/ \ bkhr / \ i / \ bkhr , 7T . . , 

hm y? T (y) = e 2 </? (0), hm ^ T (y) = e 2 <p (-), (84) 

y — >o y — ►£ 

which turns (1531 into 



?o(£ )^o(^)e- 6feftT - ^ (0)^ (0)e^ = 0. (85) 

It is clear that (1551) can be satisfied Vr iff fo(^) = = <^ (0). States with such a property 
belong to the domain of w defined by (146j) . 

In order to construct the 'evolving states' that vanish at the boundaries, consider the 
basis vectors f n {y) = e l2bny . Then, f n {y) — f m (y) satisfy the condition ( 1851) . Making use of 
flHTj) we get 

/i-e*ta(|\ 2 " / l + tan2(%)~ 

7n^,Tj ^• + e 6Mr tan (|) y l y e -bkhr _|_ g&Mr tan 2 (y) ' 1 ' 

where / n (?/,T) := fn,r(y)- Moreover we have 

(87) 



J <V 7 2 V y 2 WBW ; l + e 26fc ^tan 2 (|) V l + tan 2 (|: 

Using the substitution x = tan(^) we get 

(/m| - ia-^-fn) = 
dy 

j _ e bkTiT x \^{n-m) f-bktvr _ ^>ktir^ x 



I ■ ~ i JokhT™ ) ( c -bkhr I c bkhr^2\2 ^ X (^) 





- / i _ e 6fcfir a .N2(n-m) x + ^ 



lf ' / '' ' Q ^ + e bfcfir x J (g-bfeftr _|_ e bkhr x 2y dX ' 



Another substitution 2; = e x leads to 



(/ml - ia-^-fn) 
dy 



2(n—m) 

I — Z x 



-ia(e- 6fc?lT - e bkhT ) / ( - — - ) - — ^r-r dz (89) 

Vo V^ + ^y (i + ^ 2 ) 2 

+Aanl \—z (1 + , 2 ) 2 dZ 



Finally, we obtain 



lf\ _ ia—f n ) = { Mn-n^-A 1 ~ 8n ( n ~ m )) sinh(6fcfrr), n^m 
m <i|/ n \ ia smh(bkhr) + 27ma cosh(6/cfrr), n = m. 
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Now, let us introduce g n m(y, T ) '■= /j^ 1 ^ so that \\g n m\\ = 1- One has 

v -f 

(flwl — ia,—g nm ) = (n + m)ab cosh (bkhr) = — — — A cosh(bkhr) . (91) 
ay 2 

The expectation value of the operator f l9Tj) . defining the volume operator, is similar to 
the classical form fl29|) . The vectors g nm may be used in the construction of a basis of the 
space of states such that <Po(j) = = <£>q(0). 



IV. SUMMARY AND CONCLUSIONS 

Turning the Big Bang into the Big Bounce in the Bianchi I universe is due to the modi- 
fication of the model at the classical level by making use of the loop geometry (in complete 
analogy to the FRW case). The modification is parameterized by a continuous parameter A 
to be determined from observational astro-cosmo data. 

The reduced phase space of our system is higher dimensional with nontrivial boundaries. 
This requires introducing a few new elements, comparing to the FRW case, into our method: 
1. Our parameter does not need to coincide with the scalar field, commonly used in loop 
quantum cosmology, and it simplifies the form of a volume operator. 2. We introduce the so- 
called true Hamiltonian. It generates a flow in the family of volume quantities, enumerated 
by the evolution parameter. It proves an independence of the spectrum of the volume 
operator on the evolution. 3. We divide the phase space of the system into two distinct 
regions: Kasner-like and Kasner-unlike. 4. We identify domains, spectra and eigenvectors 
of self-adjoint directional volumes and total volume operators in the Kasner-unlike case. 5. 
We identify the peculiarity of the Kasner-like case due to complicated boundary of the phase 
space region. We propose to overcome this problem by dividing this region further into 
smaller regions, but with simpler boundaries. 6. Given a small subregion for Kasner-like 
case, we propose a canonical redefinition of phase space coordinates in such a way, that we 
can arrive at relatively simple form of volume operator and at the same time can simply 
encode the boundary of the region into the Schrodinger representation. Then, from a number 
of different operator orderings we chose the simplest one. We find domain, spectrum and 
eigenvectors of the volume operator. The spectrum is given in an implicit form in terms of 
special functions. 7. Having the true (self-adjoint) Hamiltonian, we introduce an unitary 
operator with an evolution parameter. 

The spectrum of the volume operator, parameterized by A, is bounded from below and 
discrete. An evolution of the expectation value of the volume operator is similar to the 
classical case. We have presented the evolution of only a single directional volume operator. 
One may try to generalize this procedure to the total volume operator. Analyzes, in the case 
of the Kasner-like dynamics, are complicated and will be presented elsewhere. 

Discreteness of space at the quantum level may lead to a foamy structure of spacetime 
at the semi-classical level. A possible astro-cosmo observation of dispersion of photons 
travelling over cosmological distances across the Universe might be used to determine the 
value of otherwise free parameter A. The discreteness is also specific to the FRW case jsj. 
The difference is that in the Bianchi I case the variety of possible quanta of a volume is much 
richer. On the other hand, the Bianchi type cosmology seems to be more realistic than the 
FRW case, near the cosmological singularity. Thus, an expected foamy structure of space 
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may better fit cosmological data. Various forms of discreteness of space may underly many 
approaches in fundamental physics. So its examination may be valuable. 

Quantum cosmology calculations are plagued by quantization ambiguities. In particular, 
there exists a huge freedom in ordering of elementary operators defining compound observ- 
ables, which may lead to different quantum operators. Classical commutativity of variables 
does not extend to corresponding quantum operators. Such ambiguities can be largely re- 
duced after some quantum astro-cosmo data become available. Confrontation of theoretical 
predictions against these data would enable finding realistic quantum cosmology models. 

Our nonstandard loop quantum cosmology method, successfully applied so far mainly to 
the FRW type models, seems to be highly efficient and deserves further development and 
application to sophisticated cosmological models of general relativity. 



Appendix A: Unitarily non-equivalent volume operators 



In both Kasner-like and Kasner-unlike cases, we have reduced the Hilbert space by re- 
moving the double degeneracy of eigenvalues for the volume operators (see the discussion 
after equations f|44l) and flSTl) ). We have used the 'natural' condition that the wave function 
should vanish at the boundaries of an interval. However, there are also other mathematically 
well-defined choices for the boundary conditions. We will demonstrate this non-uniqueness 
for the Kasner-unlike case. Similar reasoning applies to another case. 

Let us begin with the equation 



d 2 

2 a p 2; 
- a d^ f=UL 



ye (0,7r/&), 



which has the solution 



/„ = Nj_ sm(-y) + N 2 cos(-y), N lt N 2 e C, 
a a 



(Al) 



(A2) 



for each value of v 6 M.+ ( v i— >■ —v does not produce any new space of solutions). Our task is 
the determination of self-adjointness of w : = 



a 2 -£p and removing the double degeneracy 



of eigenvalues. The symmetricity condition reads 



If" = If 



■w/b 



f'f 



n/b 



+ / f'f- 



(A3) 



We can set: 
. /(0) = /(tt/6) = = 
. f (0) = f'(n/b) = 
. /(0) = /'(tt/6) = : 



/„ = sin(-2/), v = ab, 2ab, 3ab, . . . 

> f u = cos(-y), v = 0, ab, 2ab, 3ab, . . . 

> f v = sin(^), v 



\ab, |a6, |afe, . . . 



where ab 
spectra. 



12iihG'-fX. All these choices are non-equivalent, since they lead to different 



16 



Appendix B: Standard quantization 

Let us change the coordinates of the Kasner-like sector phase space (fii,^), defined by 
f l5T]) . into a new canonical pair as follows 

X := ^2Sl[ and P := Q 2 V^h, (Bl) 

where 

(X, P) e (0, V2di/6) x R, {X, P} = 1. (B2) 
In the new variables the volume (1531) reads 

1 Vl = ip 2 + ix 2 . (B3) 



AttG^X 2 2" 

Thus, in these variables the volume has a form of the Hamiltonian of the harmonic oscillator 
in a 'box' (0, v/2di/6). 

In the Schrodinger representation, i.e. X := x and P := —ihd x , a standard quantization 
yields 

1 ti 2 1 

— — u = + -x 2 , (B4) 

which corresponds to the 'nonstandard' quantization (160]) with the parameters m = k = 1/4 
and |/ = a/2x (with ft = 1). 

Thus, we can see that the prescription defined by ( 155|) and (1551) includes not only a 
standard prescription, but many others. We have completed only one, corresponding to the 
well known harmonic oscillator, as an illustration. 
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